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Giant Ambipolar Rashba Effect in a Semiconductor: BiTel 
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We observe a giant spin-orbit splitting in bulk and surface states of the non-centrosymmetric semiconductor 
BiTeL We show that the Fermi level can be placed in the valence or in the conduction band by controlling the 
surface termination. In both cases it intersects spin-polarized bands, in the corresponding surface depletion and 
accumulation layers. The momentum splitting of these bands is not affected by adsorbate-induced changes in 
the surface potential. These findings demonstrate that two properties crucial for enabling semiconductor-based 
spin electronics - a large, robust spin splitting and ambipolar conduction - are present in this material. 

PACS numbers: 73.20. At, 71.70.Ej,79.60.Bm, 79.60.-i 



The relativistic spin-orbit interaction (SOI) lifts the usual 
Kramers spin degeneracy in electron systems that lack inver- 
sion symmetry. It lies at the origin of many subtle and interest- 
ing effects in the electronic structure of materials such as the 
emergence of topological insulators (TI), a new quantum state 
of matter. In the bulk of materials with non-centrosymmetric 
structures, such as the zincblend and wurzite structures, it 
gives rise to the Dresselhaus 1 1 ] and Rashba Q effects. An 
analogous effect, the Rashba-Bychkov effect, describes the 
lifting of the spin degeneracy at surfaces and at asymmet- 
ric interfaces, where inversion symmetry is also broken l3l . 
The SOI is a general phenomenon, but it is especially rel- 
evant in solids containing high-Z elements because of their 
large atomic spin-orbit parameter. The characteristic split- 
ting in energy and momentum was first directly observed by 
angle-resolved photoelectron spectroscopy (ARPES) on the 
Au(lll) surface (U. The predicted polarization of the elec- 
tronic states was confirmed by spin-polarized ARPES EJO, 
and the Rashba scenario has been extended to other surfaces 
and interfaces CHIS1- 

The vision of an all-electric control of spin transport in new 
device concepts explains the strong current interest for mate- 
rials with large Rashba or Dresselhaus effects. Future devices 
operating at room temperature will require a large separation 
between the spin-polarized bands and the ability to tune the 
position of the chemical potential over a broad energy range. 
Whereas the former have been reported in surface alloys with 
high-Z elements such as Pb or Bi, only limited tunability has 
been achieved so far. 

BiTel is a non-centrosymmetric semiconductor for which 
theory predicts a large bulk Rashba effect, and the emergence 
of a topological insulating phase under pressure [16]. Ishizaka 
et ah 1 17 ] used spin-resolved ARPES to reveal spin-polarized 
states with a large momentum splitting. They assigned them 
to a quantum- well state (QWS) confined in the accumulation 
layer that appears because of band bending in the surface re- 



gion. This interpretation has been questioned in part by more 
recent ARPES data and theory that show the coexistence of 
surface and bulk bands near the Fermi level 1X51 . 

In this Letter, we show that in BiTel the chemical potential 
can be moved well into the conduction band or the valence 
band by controlling the surface termination. Remarkably, a 
giant spin splitting at the Fermi surface is observed in both 
cases. First-principles relativistic calculations indicate that 
both the surface and the bulk bands are split by the SOL We 
also prove that size of the Rashba effect is largely insensitive 
to changes in the surface potential. Therefore, the splitting 
has mainly an atomic origin. These results establish BiTel as 
a versatile material, characterized by the coexistence of very 
large ambipolar bulk and surface Rashba effects. 

We performed ARPES experiments at the Electronic Struc- 
ture Factory, beam line 7.0.1 of the Advanced Light Source. 
The energy and momentum resolution of the hemispherical 
Scienta R4000 analyzer were 30 meV and 0.1°. High quality 
single crystals of BiTel, in the form of platelets, were grown 
by chemical vapor transport and by the Bridgman technique, 
and characterized by x-ray diffraction and transport. They 
showed a metallic conduction due to a small (<2%) devia- 
tion from stoichiometry. The samples were mounted on a He 
cryostat and cleaved in UHV to expose flat, shiny surfaces. 

First-principles electronic structure calculations were per- 
formed within the density functional theory (DFT) framework 
employing the generalized gradient approximation (GGA) as 
implemented in the Quantum-ESPRESSO package |[T9l . 
Spin-orbit effects were accounted for using the fully relativis- 
tic norm-conserving pseudopotentials acting on valence elec- 
tron wavefunctions represented in the two-component spinor 
form [20]. The surface band structures were obtained using 
a slab model consisting of 39 atomic layers. Since the BiTel 
crystal has no inversion symmetry, the surfaces of the slab 
are necessarily different. The model considered includes two 
experimentally relevant terminations (Te and I). The surface 
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FIG. 1: (Color online) (a) I 4d, Te 4d and Bi 5d core level spectra measured at his — 120 eV on Te- (red, top) and I-terminated (green, 
bottom) surfaces. (b,c) Schemes of surface band bending for the two surface terminations, (d) ARPES intensity map along the TK direction 
measured at 93 eV and T = 40 K, compared to (e) the projected slab band structure calculated from first principles. The size of red symbols 
is proportional to the magnitude of projection onto the surface Te atoms indicated in the inset. The continuum of bulk states is shown in blue. 
(f,g) Corresponding plots for the I-terminated surface. The size of the green symbols is proportional to the contribution of the surface I atoms 
indicated in the inset. 



bands were extracted from the slab band structure by project- 
ing the Kohn-Sham wave functions onto atomic wave func- 
tions at the surface layer. The bulk and the slab band struc- 
tures were aligned by matching the potential in the middle of 
the slab with the bulk potential. 

BiTel has a trigonal layered structure, with Bi, Te and I 
planes alternating along the c axis. The Bi and Te planes 
are covalently bonded to form a positively charged (BiTe) + 
bilayer. The ionic coupling between the bilayer and the ad- 
jacent I - plane defines the natural cleavage plane [21 j. The 
topmost layer - Te or I - is identified by the relative intensi- 
ties of the Te and I Ad core levels, as in Fig. 1(a). Ideally, due 
to the lack of inversion symmetry, the surface termination is 
uniquely determined by the direction of the c-axis. However, 
repeated cleaves of the same crystal randomly expose both 
terminations due to the occurrence of stacking faults, which 
also explains the observation of 6-fold symmetry in the Laue 
patterns (not shown). We have measured several 'pure' sur- 
faces and on 'mixed' surfaces that presented areas with both 
terminations l22l . Data for the former are illustrated in Fig. 1. 
The surface charge - positive for Te, negative for I - induces 
band bending in opposite directions for the two terminations. 
The Fermi level lies into either the conduction or the valence 
band, giving rise to a charge accumulation or, respectively, 
depletion layer. This is schematically illustrated in Fig. l(b,c) 
and substantiated by the ARPES data. 



Figure 1(d) is an ARPES intensity map of the Te-terminated 
surface, measured along the TK high- symmetry direction 
(FK = 0.82 A" 1 ) of the Brillouin zone. The most promi- 
nent feature is the split parabolic band (SSx e i) straddling the 
Fermi level E F . The two subbands have minima at —0.32 eV 
and are offset by ±0.05 A around T. This is consistent 
with previous data, and with a Rashba interaction one order 
of magnitude stronger than for the Au(lll) benchmark case 
ifTvTl . This feature is well reproduced in the first-principles 
band structure of Fig. 1(e). It shows that the spin-split state 
is localized in the topmost bilayer, and partially overlaps with 
conduction band states that exhibit a smaller momentum off- 
set. The bulk signal is too weak to be identified in Fig. 1(d), 
but it can be discerned between SSx e i and Ep at closer in- 
spection (22). The projected bulk valence band exhibits gaps 
supporting other surface states. These states exhibit an even 
larger, but previously unnoticed, splitting. Their spin polar- 
ization has non negligible radial and out-of-plane components 
(22), at variance with the simple Rashba scenario, and simi- 
lar to recent observations on topological insulators (23). Only 
part of this manifold, labeled SSt C 2, is visible in the ARPES 
map, due to matrix elements and to the overlap with the bulk 
continuum. One can readily recognize symmetrically split 
subbands, with maxima at ±0.2 A -1 and —1.3 eV, crossing 
at T at — 1.8 eV. Bulk valence band states are visible between 
SSxei and SSx e 2, with a maximum 0.37 eV below the bottom 
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FIG. 2: (Color online) (a) ARPES intensity from a Te-terminated 
surface, measured at ^ 0.3° off normal emission, in a photon energy 
scan between 84 eV and 162 eV, plotted as a function of k±, the 
wave vector along the c-axis. Both the upper and the lower branch of 
SSxei are visible. The white dashed curve is the calculated dispersion 
of the corresponding bulk valence state, (b-d) Constant energy con- 
tours measured at the energies marked by the corresponding (b,c,d) 
horizontal lines in panel (a). 



of SSxei- These states are also offset in momentum, reflecting 
a large bulk spin-orbit splitting. 

The picture from the I-terminated surface (Fig. l(f,g)) is 
quite different. The electron pockets around T are replaced by 
hole pockets from a spin-split state (SSn) with a strong pro- 
jection on the surface I atoms. The momentum off est is again 
quite large, of the order of ±0.2 A -1 . A precise determination 
is difficult because the top of the band lies above Ep. The top 
of VB is also located above Ep. There is a complete change 
from electron to hole carriers with respect to Fig. 1(d), which 
demonstrates ambipolar conduction in BiTel. The total change 
in band bending between the Te- and I-terminated surfaces, 
estimated from core level spectra (22) . is AEbb = 0.9 eV, 
to be compared with the estimated energy gap AE g ~ 0.38 
eV ifTTl . Here again, additional surface features, partially 
overlapping the valence band continuum, can be identified at 
higher energy both in the ARPES map (SS12) and in the cal- 
culations. 

The assignment of the spectral features of Fig. 1 to sur- 
face or bulk states is further supported by the data shown in 
Fig. 2(a). It illustrates the photoelectron intensity measured at 
near-normal emission from a Te-terminated surface as a func- 
tion of fc_L, the wave vector along the c-axis. Bands SS Te i and 
SSTe2 exhibit an intensity modulation but no dispersion, as ex- 
pected for true surface states. By contrast, state VB exhibits a 
~ 0.7 eV dispersion along TA, consistent with its bulk char- 
acter, and well reproduced by the calculation (dashed line). 
This is confirmed by constant energy maps (CEM) of these 



states. The CEM measured at E F for SSjei (Fig. 2(b)) has two 
concentric contours typical of the Rashba scenario. The ex- 
ternal contour, warped by the interaction with the lattice po- 
tential, has a 6-fold symmetry, consistent with time-reversal 
symmetry for a surface state (24). In the CEM measured at 
the crossing point of SSt C 2 (Fig 2(d)), the inner contour has 
collapsed to a point at T. The outer contour again exhibits 
a 6-fold symmetry. Interestingly, it is not closed around T, 
but it is broken into 6 disconnected pockets aligned along the 
6 equivalent TM directions. This can be seen as the limit of 
strong warping, reflecting a large in-plane asymmetry of the 
surface potential (25) . Conversely, the CEM through VB (Fig. 
2(c)) shows a single 3-fold contour. This symmetry reduction 
is not due to partial extinction of a 6-fold contour induced by 
ARPES matrix-elements, because the pattern remained locked 
to the crystallographic directions when the crystal was rotated 
around the surface normal. Therefore, the 3-fold pattern of 
VB reflects the 3 -fold symmetry of the bulk potential, and 
confirms the bulk character of this state. 

We now turn our attention to the origin of the very large 
Rashba splitting. Competing models beyond the standard 
Rashba scenario have been proposed. They alternatively stress 
atomic contributions (26), the in-plane anisotropy of the sur- 
face potential [27], the asymmetry of the wavef unctions (28) , 
or the local orbital angular momentum (29). A bulk origin has 
been invoked for BiTel [17], but this conflicts with the surface 
nature of the relevant states. 

In the standard Rashba scenario the size of the splitting is 
controlled by the gradient of the surface potential, and this 
prediction was found to be consistent with the properties of 
QWS formed in an accumulation layer at the surface of the 
TI Bi2Se3 [30]. In order to test this hypothesis for BiTel we 
have changed in a controlled way the surface band bending, 
and hence the gradient of the potential in the surface region. 
This was achieved by depositing increasing amounts of potas- 
sium on a Te-terminated surface. Adsorbed K atoms donate 
electrons to the CB, leaving a positively charged surface layer 
which enhances the downward surface band bending. Figure 
3(a) illustrates the evolution of SSx e i as a function of K cov- 
erage. Movies of complete K dosing experiments for both Te- 
and I-terminated surfaces are available in Ref . [22j 

As expected, SSTei » as well as all core levels (22), shift to 
lower energies, following the change in band bending. The 
total shift at saturation K coverage is 0.12 eV, bringing the 
bottom of SSTei band 0.44 eV below Ep. A closer inspec- 
tion shows that the energy shift of SS Te i is rigid. Figure 3(b) 
shows the Fermi wave vector &p of the outer branch and the 
momentum offset ko. They were estimated from a parabolic 
fit of the dispersion, keeping the effective mass unchanged. 
Within error bars ko remains constant. The bottom of SSTei 
gives an upper limit for band bending. Therefore a >40% 
change in the surface band bending has no measurable effect 
on the strength of the Rashba effect. This experimental ob- 
servation strongly suggests that other parameters, namely the 
atomic spin-orbit parameter of the heavy elements, determine 
the large spin splitting. Figure 3(c) shows that the surface car- 
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FIG. 3: (Color online) (a) Evolution of the spin-split SSiei surface 
state band as a function of K coverage. The red parabolas are the 
result of a fit to the band dispersion along TK. (b) The Fermi wave 
vector fcp of the outer branch (solid symbols) and the Rashba splitting 
ko (empty symbols) are plotted as a function of the measured surface 
band bending, (c) The corresponding surface electron density. 



rier density, estimated from the area of the electron pockets, 
varies linearly with the downward shift of the split bands. This 
is again consistent with a constant fen- By contrast, deviations 
from linearity have been observed for the QWS at the Bi 2 Se3 
surface t30l . 

In summary, we have shown that large ambipolar bulk and 
surface Rashba effects coexist in the non-centrosymmetric 
semiconductor BiTel. The Fermi level at the surface lies ei- 
ther into the valence or the conduction band, depending on 
the nature of the topmost layer. The surface termination is ran- 
domly chosen by cleavage due to stacking faults in bulk crys- 
tals, but it could be precisely determined in artificially grown 
thin films. It is realistic to consider that regions with opposite 
band bending could then be patterned on a substrate, opening 
new perspectives for the manipulation of spin-polarized states. 
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